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GROUND-WAVE PROPAGATION AT MEDIUM FREQUENCY IN BUILT-UP AREAS 
J.H. Causebrook, Ph.D., B.Sc, C.Eng., M.I.E.E. 



1. Introduction 

For more than 50 years concern has been expressed 
about the loss of field strength experienced in built-up 



areas aL m.i. 



I lie concerf 



iiiwicdocLj idLdty 



U.K. as a result of the poor fields received from local broad- 
cast transmitters in cities, from transmitter sites on the out- 
skirts. A study has therefore been made with the objectives 
of understanding the nature of this loss, and producing a 
method of predicting its extent. 

It is well known that higher frequencies in the m.f. 
band have smaller ground-wave field strengths for a given 
transmitter radiated power, but it was suspected that in 
towns the difference was greater than conventional theories 
would predict. Many m.f. field strength measurements 
have been taken over the years but they were not considered 
adequate for the present analysis. This was because they 
were mostly taken to determine the extent of service areas 
and did not provide a continuous series of field measure- 
ments from close to the transmitter to beyond the town to 
be served. 

Thus the first part of the project consisted of a 
measurement campaign. The results confirmed that 

losses at higher frequencies exceeded the values expected 
from existing theory and further inconsistencies were 
found. The second part of the work consisted of a 
theoretical study, and an extension of the theory of ground- 
wave propagation was developed which could cover the 
practical situation of built-up environments. 

The general form of the extended theory was written 
as a computer program with certain weighting factors left 
undefined. These factors were adjusted by an optimization 
procedure within the program so that a satisfactory fit to 
the measurements was achieved. The resulting values for 
the weighting factors indicated some imperfections in the 
theory. A method of overcoming these is suggested in this 
report, but the detail of this is left for future work. 

The final part of the project was to write another 
computer program which can be used for prediction work. 
This program not only uses the new theory to allow for the 
built-up environment but also an integral equation method 
of calculation which is more satisfactory for profiles which 
are irregular and inhomogeneous. The comparisons be- 
tween the results of this program and the measurements 
taken are good. 



2. The measurement campaign 

The measurement campaign was concentrated along 
two radial lines from the Brookmans Park transmitter 
(29 km to the north of central London), both passing 
through the London area. These were on bearings of 181 



and 169 degrees east of grid North, extending to about 50 
km from the transmitter. Routes for the measuring vehicle 
were chosen which ran as close to the radial lines as 
possible. The measurements were made on three fre- 

„.,„.^»;nr- onQ iniA ~.r,M iaci i^u, +Ucc:n i^^in^ +v,q Ban 
^ut^i luicS . ijKJOf 1^1^ anu i-tD/ txi iZ, li icoc i^cuig (.1 ic DDvy 

programmes Radio 4, Radio 1 and Radio London respec- 
tively. 

The measurements were made with equipment which 
was based upon that described in Reference 1. However, 
modifications and additions were incorporated to produce 
a more satisfactory arrangement for the particular task. 
Firstly, a pen recorder was incorporated so that a record 
could be obtained which had a chart displacement propor- 
tional to distance travelled (using pulses from the vehicle 
speedometer). Secondly, the original receiver was unsatis- 
factory because of non-linearity and instability of gain and 
frequency response, so it was replaced by a better instru- 
ment. Thirdly, the original equipment was meant to take 
readings while the vehicle was stationary, so it was necessary 
to overcome the directional property of the loop aerial by 
manually keeping it on bearing. 

The chart recordings were highly variable, mainly due 
to multipath propagation. To obtain a more convenient 
set of results the charts were analysed in sections, equiva- 
lent to 500 metres along the radial line. From this 
analysis median values of magnetic field strength, H, were 
obtained. It is usual to quote results in electric field 
strengths, E, but the conversion in this case would depend 
upon an unknown space impedance. It was therefore 
preferred to retain the signal strength measurement in 
terms of // in amperes per metre. The analysed values were 
plotted on graphs of 20 \og^^(Hd)^p^ against distance, 
d^^. These graphs are reproduced in Figs. 1 and 2. 

The interesting feature of these plots is that, in the 
first section, the fields are higher than expected and then a 
steep drop occurs to be followed by an increase beyond 
the densely built-up area. An extra feature should be noted 
on the 169° radial, which is a further reduction associated 
with the tall buildings in Croydon. It will also be noted 
that the steep drop does not occur at the same distance for 
each frequency. For 908 kHz it is approximately at 30 km 
from the transmitter and for 1457 kHz it is at 20 km; the 
distance of this reduction-for 1214 kHz lies between these 
values. 

In addition to the van measurements some values 
were obtained away from the vehicle using the internal 
aerial and calibration system in the receiver. These were 
taken at distances from the transmitter of between 400 
and 2500 metres. They are plotted on Fig. 3. The purpose 
of these measurements was to obtain results which were 
not restricted to roads. Thus they could be taken in areas 
as free of man made objects as possible and the radial line 
could be followed more closely. In this way a determin- 
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Ffg. 1 - Road measurements made along the radial 18J° E.G.N, f^ree frequencies} 
• 908 kHz • 1214 kHz • 1457 kHz 



ation of the point where the graphs of 20 \Qg^^{Hd) should 
intersect the zero-distance axis will be known to higher 
accuracy. Thus in later graphs, Figs. 8 to 13, we are able 
to match calculations and measurements. 



3. Theory 

The existing calculation methods depend upon the 
assignment of values to ground conductivity, a, and, less 
significantly, relative permittivity, e. Then the Sommerfeld 
theory gives a method for prediction of field strength over 



a homogeneous smooth earth. A description of such a 
calculation is given in Reference 2, and from this method a 
set of curves can be drawn for simple predictions. If 
there is a limited change of ground constants along the path 
then a prediction could be made by the Millington method. 

For towns it has been the practice to assign an 
effective conductivity value. This procedure was tried in 
an attempt to reproduce the measurements mentioned in 
the previous section, but the results were poor. Thus the 
possibility of producing more satisfactory methods was 
explored. " " 
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Fig. 2- Road measurements made along the radial 169 E.G.N, (tiiree frequencies) 
• 908 kHz •1214 kHz • 1457 kHz 



Hufford gives an integral equation which permits 
calculation of field strength over smoothly varying inhomo- 
geneous terrain. The ground-wave attenuation factor g{/?) 
at distance R relative to a perfectly conducting plane earth 
is determined in this method by the equation: 



gt^j=l 



- /- (t? + il') e-'^f g(r} / dr (1) 



2-n 
where^ = — <X = w/ave length) 
X 

<^ = angle (in radians) between the line from the 

receiver to the surface at range r and the tangent 

to the surface, at this point. 
5 = distance between terminals* subtracted from the 

sum of distances from terminals to the surface at 

r (see Fig. 5, Section 4). 
1? = relative surface impedance atr, 

•The ward terminals refers to the transmitter and the receiver at 
point J?. 
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Fig. 3 - Measurements away from the vehide 



For an uncluttered surface and a plane wave incident at an 
angle Q to the normal to the surface we have: 



[e-i60aX-sin^ei'''" 
6 - j6DoX 



(2J 



In a built-up area it will be as if another medium were 
placed between the ground and the 'free-space' above. 
Thus we require a modified value of t? appropriate to these 
conditions. If we were to think of this modified value as 
resulting from an 'effective' conductivity we would be 
restricted to a value of 17 which has roughly equal real and 
imaginary parts. An initial attempt to match the measure- 
ments this way indicated that what was required was a value 
of t? which had a much larger positive imaginary part. The 
reasons for this can be appreciated from an examination of 
the Appendix (Section 10). 

The two sections that follow give a theoretical 
plausibility to the proposal to adopt a modified value of 17 
as suggested by the empirical work. The analyses give an 
indication as to the form and magnitude of 7? that is appro- 
priate to built-up areas. Exact values cannot be expected 
from these analyses asr 

1, The building layer is treated as if it can be described 
by a macroscopic theory, but in practice the sizes of 
structures are not insignificant compared to a wave- 
length. 

2. Only a single plane wave is considered. 





Radial bearing 


Freq. kHz 


181° 


169° 


908 
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1214 


• 


a 


1457 


• 


B 



3. The method does not consider signals arriving by 
scattering (multipath} from points other than those on 
the line between transmitter and receiver. 

3.1. Low building density theory* 

At l.f. and m.f. many man made structures (lamp 
posts, steel framing, house-wiring and plumbing) may be 
considered as earthed parasitic unipoles. Thus an incident 
electric field induces a current in these structures. To 
analyse the effect of this we will use a co-ordinate system 
and geometry as shown in Fig. 4. It is possible to assume 
a current density in the vertical direction given by: 



/(2l = j-Jt(z)£, 



(3) 



where tj^ is the intrinsic impedance of free space and K\z) 
is a value decreasing with height and dependent on 
the length of the unipole and its characteristic 
impedance. 

From Maxwell's equations in the building layer: 



- - i-l] +K{z}] £, 



(4) 



If a plane wave exists in the space above this layer the com- 
ponent of the propagation coefficient in the >> direction 

* The original idea for this approach was suggested by P. Knight. 
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the model 
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(see Fig. 4) is -'§ %\n9 in all the media (free space, building 
layer and the ground). 

Now let us consider a small stratum in the building 
layer, between the heights z = t and z = t + bt,\n which we 
consider the vertical propagation coefficient as constant at 
7(?) denoted as F and K(z) constant having the value K(t) 
denoted as K . At t there will be downgoing and upgoing 
waves characterized by H^ and H^^. Thus in this stratum: 



//v 



H^e 



r(z-r) 



+ H..e 



-r{z-t) 



^-iPy sine jg| 



The corresponding electric field components will be: 



iiS dz iP 



H^e 



r(z-f) 



H..e 



-r(z-f) 



J7?o dH^ ??oSin0 



' pa +K') dy -i+K' 
Also from Maxwell's equations: 



//. 



-JPy sine 
(6) 

(7) 



Therefore from Equations (5) and (6) at z = t: 



Vit) 






and at z = t + 8t: 



nit + bt) 



T H^e^" - K.-^" 



rst . 



j/Ji/^e' +^ue 



^-rst 



(11) 



(12) 



Approximating for e*"^^^ and from Equations (11) and 
(12) we get: 



Vit + St) 



Vit) +— St 



1+jr?{/'))35f 



(13) 



Because both rjit) and pSt are small this may be further 
approximated to give: 



H. 



dE^ dE\, 



sin 



PVo \ 9.y dz [l +K' jS^ 



H^ 



Therefore: 



r = ji3 1 



sin 



1 +K 



(8) 



(9) 



Subject to certain criteria it is valid to accept the definition 
of relative surface impedance as being given by: 






(10) 
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r}(t + bt]=r}(t) +jj3 1 



sin 



1 +K 



\bt 



(14) 



where we have also substituted for F using Equation (9). 

If we consider the depth of the stratum reducing to an 
infinitesimal and we integrate from ground level to the top 
of the building layer we get: 



dr? = j|3 



1 +K{t) 



\6t 



(15) 



T7{0) 



J 

o 



5- 



Therefore 




(16) 



in the present problem we wish to consider waves at near 
grazing angle. To simply put_ to 90° throughout the 
above analysis would introduce some awkward zeros. Thus 
9 must be considered as only approaching 90 . 

To solve the integration in Equation (16) it is neces- 
sary to assume the form that K should take. There are 
many factors which contribute to a discussion about the 
best form for this variable, such as the dimensions of the 
re-radiating objects and the fact that a given area will have 
a mixture of these. For the present purposes a triangular 
height distribution is a reasonable compromise. Thus: 



Kit) = 



nh'-T]^ 



2Vc 



t 



(17) 



Where n = the number of structures per unit area and 

i7j, = the characteristic impedance of the unipoles. 

The predominant objects in the less built-up areas are two- 
storey houses and lamp-post like structures. They both 
have similar heights but there are more houses in most 
areas. With the use of Reference 6 we get t?^, approximately 
equal to 20 ohms for houses, and 200 ohms for lamp-posts. 
Thus the houses will dominate the value of K. It is now 
possible to assume that houses are roughly cubic in shape 
so that n/z^ equals the fraction of area covered by houses, 
B say. Thus Equation (16) becomes: 



nih) = n(o) + '§h 



1 - 



log„(1 +105) 



10fi 



(18) 



3.2. High building density theory 

The analysis of Section 3.1 becomes progressively 
less valid as the building density becomes large. This is for 
several reasons, but the most important is that the impe- 
dance of the tops of buildings is neglected. To determine 
the type of formulation required in the highly built-up 
areas consider grooves cut into a highly conductive surface 
and exposing poorer conductivity at the base. In this case 
we may distinguish areas giving three different surface 
impedance: 

1. Roof tops where we will assume i?(/z) equals zero 

2. Streets parallel to the line of propagation where t7(/z) 
equals 77(0) 



7?(/z) = 



77(0) -I- jtanjS/; 
1 + J77(o)tani3/z 



(19) 



Provided ph<^lA, only a small error is introduced if 
we make a similar approximation to that which gave 
Equation (14). 



Then 



r^m = rt(o) + \&h 



(20) 



It is not practicable to use these values individually in 
the calculation system. A compromise surface impedance 
may be obtained by taking a weighted average of roof tops 
and the two directions of streets. Thus as 5 is the fraction 
of area covered by buildings we may assume that each type 
of street covers an area Vzd - B), which leads to: 



n(h)= (1 -3)71(0) + 



1 -B 



m 



(21) 



3.3. A combined result 

It is now convenient to combine the results given by 
Equations (18) and (21). However, before doing this it 
should be pointed out that 17(0) will not necessarily remain 
at the value which would have existed before the area was 
built-up. Let us refer to this indiginous surface impedance 
as T?j. It is then possible to denote the combined result 
thus: 



T? = fj(5)T?i+f2(5)j/3/z 



(22) 



where fj (B) and f^ (B) are functions of building density. 
The precise form of these functions is not obtainable from 
theory because building density is not a perfect parameter 
and the theories are approximate. To overcome this we 
may determine the functions via the measurements. This 
may be done with the aid of a computer optimization 
process, but first it is necessary to decide on the general 
form of the functions and the constraints which should be 
placed upon these forms. For example, a non-monotonic 
situation may be expected because Equation (18) indicates 
the imaginary part increasing with building density, whereas 
Equation (21) shows a decrease. An obvious requirement 
is that the result must tend to 77, when buildings are absent. 
These considerations led to the following: 



f,(B) 
f^iB) 



(1 -B)''i 

x,VS + x,5 + x^5^-^ 



(23) 



where the Xj values are to be determined by optimization. 

4. The computer calculation system 

Before a computer program could be written it was 
necessary to decide how best to solve an integral of the form 
given in Equation (1).* Of course, data will be discrete in 
form, i.e. a quantum of data representing a distance D 
within a profile. Let the quantum be defined thus: 



3. Streets at right-angles to the propagation direction. 
Here r7(ft) may be obtained by analogy to a trans- 
mission line. From Jordan, Chapter 8, we obtain: 



The method of solving the integral equation used here is a further 
development of a method initially proposed by G.D. Monteath. 
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Q„={v„ + ^n)e-'^^^" (24) 

where the quantities are defined in Fig. 5. 



The essential part of Equation (1) which requires 
solution now becomes: 



To obtain the value Ij, when N is greater than two, Equa- 
tion (26) may be expanded as a power series. To obtain 
satisfactory accuracy in the most critical cases it was found 
that the series needed to contain the first six terms: 



N 



E 






N 



Qn 



g{r) dr 



y/r(ND-r) 

n = 1 (n-1)Z) -n 



^-=^ 

E'^"'" 



(25) 



In most intervals we can assume that g varies linearly with r. 
However, considering the immediate vicinity of the trans- 
mitter, this assumption is invalid in the first interval and is 
inadequate in the second when there are only two intervals. 
This limitation is overcome by the fact that the first 
interval should obey the Sommerfeld/Norton formula 



4{crf 
g = 1 - j(7rcr)'''' - 2cr + ]{nc^r^ f'^ + Hnc^r^ )'''= 

(29) 

where c is a constant within the interval. A value of I^ 
may now be produced by analytical means. The assump- 
tion of linearity will now give us the remaining integral as: 



g = 1 - j \/TTpe " erfc(j\/p) 



(26) 



where A„ and R^ may be derived analytically. 
Then if: 



W = 



(30) 



Thus to calculate the value, Gj, of g at the end of the first 
interval: 




(31) 



we get 



Dtt . 



To calculate G^ it is reasonable to assume: 



(27) 



W 



2DiT 



N - 1 
n = 1 



QnIn-GNANGN-1 



(28) 



(32) 



W + Qt^R^ 
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F/'g. 7 - Output of computer program 
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5. The computer optimization process 

The field calculations can be made by the methods of 
the previous section. However, the value of r}^ to be used 
still requires satisfactory numbers for the quantities x^ to 
x^ of Equation (23). This was achieved with the aid of a 
computer optimization technique based on a method by 
Dickenson^ which minimized the sum of the squared 
difference between calculation and measurement. 

The limitations on computer storage space resulted in 
the decision to optimize in six different sets; one for each 
frequency on each radial. The values obtained were com- 
patible for a given frequency, but not so between fre- 
quencies. The frequency differences varied smoothly 
across the band but were not constant as expected from the 
theoretical considerations. The differences were suf- 
ficiently large and significant to make a choice of compro- 
mise values for all frequencies in the 900 to 1500 kHz 
range impossible. The frequency difference obtained for 
the fj(5) (Equation (23)) value could be explained by a 
'skin depth' phenomenon, i.e. the 'effective' conductivity is 
different at different frequencies as a result of differing 
penetration depths into an earth, whose conductivity varies 
with depth. A similar explanation may be found for the 
fj (B) difference but this is less clear at present. 

A further possible reason for the discrepancy can be 
seen by examining Fig. 14 of the Appendix, (Section 10) 
which discusses the phenomenon, expected from theory, of 
a very low field strength at a critical range and for surface 
impedance close to a critical value. The sharply defined 
area with very low field will be less marked in practice 
due to fields being scattered from points other than those 
on the line between transmitter and receiver. In the 
optimization process a calculated field from such an area 
would be very different from the measurements and the 
optimization parameters would adjust to force the values 
away from this situation. 

Except for the above problem the agreement between 
empricism and theory was satisfactory. The optimization 
gives f j and f^ a"^ functions of frequency: 



f,(S) = (l-5)^^/^ 



f, (B) = 




- V23B + 0-35B 



1-5 



(33) 



6. A working program 

A computer program has been written* so that field 
strengths may be calculated for an operational purpose, i.e. 
to predict a service area. A flow chart of this program is 
given in Fig. 6 and an example of the output is shown in 
Fig. 7. This output gives the wavelength in metres, 
labelled as 'LAIVIBDA(IVl)'. The columns are as follows: 

1. The distance of field point from transmitter in metres. 



■ Much of the programming work was done by B. Tait. 



2. The terrain height in metres. 

3. The average height assumed for buildings in metres. 

4. The percentage of area covered by buildings. 

5. The ground conductivity in miliisiemens per metre. 

6. and 7. Q^ as defined by Equation (24). 

8. The magnitude of the loss relative to a perfectly con- 
ducting flat surface. 

9. The phase of the loss factor. 

10. Column 8 expressed as dB. 

11. The electric field strength in millivolts per metre that 
is equivalent to the calculated magnetic field in a 
plane wave in free space. 

This program was run to produce the fields which 
were described in Section 2. The results of these calcula- 
tions are shown in Figs. 8 to 13. Each figure has the 
measurements plotted on it for comparison. It will be 
seen that the calculations closely follow the measurement 
features. The calculations for 908 kHz, Figs. 8 and 11, 
repeat the feature of high fields up to approximately 24 km 
from the transmitter. The steep decrease in field is well 
matched on each diagram as is the dip near the centre of 
London. This being the case even though these features 
differ dependent on frequency (compare Figs. 8, 9 and 10). 
Figs. 12 and 13 clearly show the repeated dip for the high 
buildings in the Croydon area. The standard deviations of 
the differences between calculations and measurements are 
as follows; 

1. 1 -9 dB for 908 kHz 

2. 3 dB for 1214 kHz 

3. 3-2 dB for 1457 kHz 



7. Future work 

Sufficient work has been done to demonstrate that the 
new calculation method is on the right lines. However, 
there is much more that could be done to fill in the finer 
details. Three points in particular require further investi- 
gation: 

1. The ratio of the electric to magnetic field in the 
built-up area requires investigation. The theory 
indicates that the mean value of the electric field will 
be lower, relative to the magnetic field, than in rural 
areas. 

2. The height gain to be expected requires confirmation 
by measurement. 

3. The frequency dependence of the factors requires 
confirmation, in view of the discrepancy between the 
measurement and the theory offered. 
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There are certain features of the calculation method 
which could be improved. In particular a method could be 
devised to allow for off-axis scattering. This could simply 
take the form of a background field which would be a limit 
to how low the field would reduce in the difficult region 
shown in Fig. 14 and considered in the Appendix (Section 
10). 

The calculation system could be turned into a yet 
more economic tool for everyday planning. The data 
required can be extracted from banks such as those held for 
UHF planning.^ The building data can be approximately 
inferred from the population data banks. 

It is probable that the methods employed here could 
be extended to other frequency bands, both higher and 
lower. 



8. Conclusions 

The problem of loss of medium-wave field strengths in 
towns has been considered for over 50 years, but no ade- 
quate prediction method has been developed. This report 
provides such a method, although further work is necessary 
to put it on a firmer basis and to turn it into a planning 
tool. 
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10. Appendix 



The new calculation system proposed in this report is 
in a domain of the ground-wave phenomenon which is 
unfamiliar. It is therefore helpful to have a picture of the 
features to be expected in this domain. This led to the 
production of Fig. 14, where a set of contours is given of 
loss (dB) and phase (degrees) relative to free-space con- 
ditions for a wave travelling over homogeneous ground. 
These are drawn on the complex plane of p, the parameter 
defined by Equation (27) and on Fig. 14. This is a useful 
choice of parameter because profile lines are radials on this 
diagram, which can be seen from the definition of p. A 
practical profile cannot, in general, be determined from 
this diagram because it does not represent the situation 
where the surface impedance varies along the path. 

The lower half of the figure represents the conven- 
tional domain and it will be seen that it is very regular and 
that the contours in it all represent loss relative to free 
space (negative dB's). The upper half represents the 
situation with buildings on the profile line and it will be 



seen that it contains a sharp irregularity and that radial 
lines can cross contours of enhanced field (positive dB's). 
The radial line making an angle of about 51° with the real 
axis passes through the deep zero singularity. Thus complex 
values of the surface impedance having angles of about 
70-5° are critical. The calculations for this diagram were 
made with the formula of Equation (26) because this per- 
mits a simpler computation for the homogeneous profile. 
This makes the assumption that the mathematical form of 
Equation (26) has a physical reality in the positive quadrant 
of the parameter p. However, a check was made by 
obtaining some results from Equation (1). The agreement 
was very close; this gives reasonablej'ustification for each 
of these equations because they have largely independent 
derivations. In addition, the extension into this new 
region produces field values that match the measured 
results thus producing a further justification. 

A further understanding of Fig. 14 may be aided with 
an hypothetical example: let us take the case where un- 
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cluttered ground would give an 7? = 0-05 + jO-05 and assume 
it is evenly built-up along the whole profile, adding jO-1 to 
give: 



Then: 



dl\= 12-7|p| 



(35) 



7? = 0-05 + jO-15 



(34) 



This results in a phase of 53-13 degrees for p which 
can be represented by the dotted black line on Fig. 14. 



The example results in the attenuation factor varying as 
shown in Fig. 15. This is very similar to the measured 
field variations (of Section 2 in the body of the report), 
even though the varying density of clutter is not included. 
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